Adiabatic quantum computation and quantum phase transitions 
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We analyze the ground state entanglement in a quantum adiabatic evolution algorithm designed 
to solve the NP-complete Exact Cover problem. The entropy of entanglement seems to obey linear 
and universal scaling at the point where the energy gap becomes small, suggesting that the system 
passes near a quantum phase transition. Such a large scaling of entanglement suggests that the 
effective connectivity of the system diverges as the number of qubits goes to infinity and that this 
algorithm cannot be efficiently simulated by classical means. On the other hand, entanglement in 
Grover's algorithm is bounded by a constant. 
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Quantum adiabatic computation |lj inherently brings 
the quantum system near to a point where a quantum 
phase transition |2| takes place. Entanglement is then 
expected to pervade the system as we shall numerically 
give evidence when quantum adiabatic computation is 
applied to the Exact Cover problem, classically classified 
as NP-complete. Furthermore, scaling of the entropy of 
entanglement appears to be linear for spin systems with 
non-local interactions, sub- linear for the XX model, log- 
arithmic for the quantum Ising model 0, 0] and bounded 
for the Grover's adiabatic algorithm. In order to sub- 
stantiate these statements we first need to recall three 
established results. 

Large pure-state entanglement of the quantum regis- 
ter is a key element for exponential speed-up of quan- 
tum computation. This result has been made quantita- 
tive by Vidal who has proven that a quantum register, 
such that the maximum Schmidt number of any bipar- 
tition is bounded, can be simulated efficiently by clas- 
sical means. The measure of entanglement proposed in 
ref. [j| is E = log 2 x, where x 1S the maximum Schmidt 
number of any bi-partitioning of the state. It can be 
further proved that E > S(p), where the von-Neumann 
entropy S refers to the reduced density matrix of any 
of the two sub-partitions. If a n-qubit quantum register 
only uses little entanglement all along the computation, 
that is x = P°ly( n ) a t most, both the quantum state as 
well as the action of the quantum gates on it can be ef- 
ficiently simulated by classical means. This implies that 
exponential speed up is only possible if entanglement per- 
vades the quantum register at some point along the com- 
putation, that is, if x ~ exp(n a ), with a being a positive 
constant, which is naturally satisfied if the entropy obeys 
S(p) <~ n b , for some positive constant b. Any algorithm 
designed to exponentially accelerate a classical computa- 
tion must create exponentially large x- An exponentially 
big x is therefore a necessary, though not sufficient, con- 
dition for quantum exponential speed-up. 

On a second separate development, entanglement for 
the ground state of many quantum spin chain systems 



has been proven to scale at quantum phase transitions 
HQ (see also [10). 

The entropy associated to tracing 
out all but L spins out of an infinite spin chain displays 
logarithmic scaling controlled by the central charges, c 
and c, classifying the universality class of the phase tran- 
sition la: 



c c + c. T 
Sl = —z- log 2 L . 
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(1) 



Results from field theory suggest that d-dimensional spin 
systems should display a leading scaling behavior com- 
pletely determined by the area of the region separating 
the partitioning of the system. For instance, when sepa- 
rating the system in the interior and exterior of a sphere 
of radius R and assuming an ultraviolet cutoff xq, the 
entropy of e.g. the interior is 



d-l 



(2) 



where c\ corresponds to a known heat-kernel coefficient 
[j| . This leading scaling behavior can be cast in terms of 
the number of spins in the system as 



S 



n d 



(3) 



Entanglement only saturates for non-critical quantum 
systems in one dimension 0. 

The third element we need to introduce corresponds 
to the quantum adiabatic computation framework intro- 
duced by Farhi et al. Q. The quantum register is ini- 
tially prepared on the ground state of a known initial 
Hamiltonian Hq . The system is then made to evolve adi- 
abatically from this Hamiltonian to a new one Hp whose 
ground state codifies the solution to an e.g. NP-complete 
problem 



H(s(t)) = (1 - s(t))H + s(t)H f 



(4) 



Slow evolution from s(t = 0) = to s(t = T) = 1 guar- 
antees that the system will not jump from the instan- 
taneous ground state of the system to the first excited 
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state. Quantum adiabatic computation is efficient pro- 
vided that the minimum gap along the adiabatic evolu- 
tion is only polynomially small in the number of qubits. 

It follows from the above arguments that quantum adi- 
abatic computation can be viewed as a time evolution in 
which there is a flow along the parameter space defin- 
ing the Hamiltonian. At a given point s c the Hamilto- 
nian approaches a quantum phase transition, character- 
ized by a vanishing energy gap. Exponential speed-up 
needs large entanglement which is also expected at some 
quantum phase transitions as discussed previously. A 
quantum computer programmed to find the solution to a 
given problem using adiabatic evolution does in fact cor- 
respond to a system that passes near a quantum phase 
transition. A quantum computation is thus equivalent to 
the simulation of a very specific quantum phase transi- 
tion. Reversely, simulating a quantum phase transition is 
known to be in general a hard problem that is in principle 
efficiently solved by adiabatic evolution if the energy gap 
does not vanish exponentially with the number of qubits. 

We shall give support to the above picture by ana- 
lyzing the span of entanglement along a quantum adia- 
batic computation applied to the Exact Cover problem, 
closely related to the 3-SAT NP-complete problem. We 
shall indeed see that entanglement seems to span over ex- 
ponentially many states in the computational basis and 
therefore the algorithm may be hard to simulate in an ef- 
ficient way using a classical computer. This is a necessary 
(thought not sufficient) condition for quantum exponen- 
tial speed-up which is apparently successfully verified in 
our case. 

The NP-complete Exact Cover problem is a particu- 
lar case of the 3-SAT problem and is defined as follows: 
given the n boolean variables {%i}i=i,...m %i = 0, 1 V i, 
where i is the bit index, we define a clause C involving 
the three bits i, j and k by the constraint Xi+Xj+x^ = I. 
There are only three assignments of the set of variables 
{xi, Xj ;, Xk} that satisfy this equation, namely, {1,0,0}, 
{0, 1,0} and {0,0, 1}. An instance of Exact Cover is a 
collection of clauses which involves different groups of 
three qubits. The problem is to find a string of bits 
{xi, X2 ■ ■ ■ , x n } which satisfies all the clauses. 

The Exact Cover problem can be mapped to finding 
the ground state of a Hamiltonian Hp in the following 
way |9j: given a clause C define the Hamiltonian associ- 
ated to this clause as 

ffc = g((l + <)(l + ^)(l + 4) 
+ (I - af)(l - cxf )(1 - ai) + (I of )(1 - ^)(1 + o%) 
+ (1 - af)(l + a|)(l - a%) + (I + of )(1 - erf )(1 - a%)) 

(5) 

where cr z |0) = |0), a z \l) = — 11). The quantum states of 
the computational basis that are eigenstates of He with 
zero eigenvalue (ground states) are the ones that corre- 



spond to the bit string which satisfies C, whereas the rest 
of the computational states are penalized with an energy 
equal to one. The problem Hamiltonian is constructed 
as the sum of all the Hamiltonians corresponding to all 
the clauses in the instance, 

C G instance 

The ground state of this Hamiltonian corresponds to the 
quantum state whose bit string satisfies all the clauses. 
The original problem stated in terms of boolean logic has 
been cast into the hard task of finding the ground state 
of a spin system with non-local two and three body inter- 
actions. The couplings depend on the particular chosen 
instance and, therefore, the spin system has not an a 
priori well defined dimensionality neither a well defined 
lattice topology, in contrast with the usual spin models 
(e.g. the anti-ferromagnetic Heisenberg model in a 2- 
dimensional square grid). This intrinsic relation between 
physical systems and difficult computational problems is 
a well established fact. For example, the ground-state 
search for some spin Hamiltonians defined on regular 
two-dimensional planar cubic lattices is proved to be an 
NP-hard problem Notice that the systems consid- 

ered in this paper, namely those arising from the Exact 
Cover problem, differ a lot from those of as l° n g 
as the detailed structure of the Hamiltonian is consid- 
ered. Nevertheless, it is a remarkable fact that the Exact 
Cover problem can be mapped to a ground-state search 
of a spin-system, as in , which provides some physical 
intuition. 

Adiabatic evolution is carried by the s-dependent 
Hamiltonian H(s) as a linear interpolation between an 
initial Hamiltonian Ho and Hp: H(s) — (1 — s)Hq + sHp, 
where the initial Hamiltonian Hq can be taken as a mag- 
netic field in the x direction 

ffo = E|(l-^). (7) 

i=l 

where dj is the number of clauses in which qubit i ap- 
pears. The ground state of H a is an equal superposition 
of all the possible computational states. Observe that 
H(s) is, apart from a constant factor, a sum of terms 
involving local magnetic fields in the x and z direction, 
together with two and three-body interaction coupling 
terms in the z component. 

Our numerical analysis is based on the random gener- 
ation of 300 instances for Exact Cover with only one pos- 
sible satisfying assignment for n = 6 up to n — 20 qubits. 
We produce the instances by adding clauses at random 
until there is exactly one satisfying assignment, starting 
over if we end up with no satisfying assignments. Ac- 
cording to , these are believed to be difficult instances 
for the adiabatic algorithm. For every instance, we have 
constructed its corresponding interpolating Hamiltonian 
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and found the ground state for s — to s — 1 in steps of 
0.01. We then consider a repartition of the system into 
two blocks of n/2 qubits and calculate the entanglement 
entropy between the two blocks as a function of s. We 
have explicitly checked on some instances that all pos- 
sible partitions produce entropies of the same order of 
magnitude (as expected from the non-locality of the in- 
teractions) and chosen to work with the first n/2 versus 
the rest. 

The results we find for the scaling of entanglement 
seem to agree with the idea that the system approaches a 
quantum phase transition along its adiabatic evolution. 
For each of the randomly generated Hamiltonians we ob- 
serve a peak in the entanglement entropy around a criti- 
cal value of the parameter s c ~ 0.7. The average entropy 
shape over the 300 instances is represented in Fig. 




Figure 1: Average over 300 instances of the entanglement 
entropy between two blocks of size n/2 as a function of the 
parameter s controlling the adiabatic evolution. A peak ap- 
pears for s c ~ 0.7. The plot also shows the increase of the 
peak as the number of qubits grows n = 10, 12, 14. 

In order to analyze how entanglement scales at the 
critical value s c , we plot the maximum entropy of entan- 
glement as a function of the number of qubits, both for 
the worst case and for the average over the 300 instances. 

The numerical analysis shown in Fig. [3 apparently 
agrees with linear scaling and matches the expectation 
that the Exact Cover problem can be viewed as a spin 
system with highly non-local three-body couplings, and 
therefore high effective dimensionality. The points can 
be fitted by a function of the type E(n) ~ .1 n. This be- 
havior would correspond to a nearest neighbor-like cou- 
pling in d ~ n dimensions, thus diverging as n goes to 
infinity. We note at this point that the evidence of large 
entanglement present in the ground-state of the system 
does not say anything about the efficient performance of 
the quantum adiabatic algorithm. Despite involving a 
highly quantum-correlated system, the running time of 
the algorithm would still be inefficient if the gap were 
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Figure 2: Scaling of the entanglement entropy both for the 
worst case and for the average over 300 randomly generated 
instances. Error bars give 95 per cent confidence intervals. 
The behavior appears to be linear. 



exponentially small. This does not seem to be the case, 
as we shall see, according to our simulations. 

Our numerical analysis is also consistent with the work 
of Farhi et al. in [|J where the minimum energy gap ap- 
pears to decrease as g m in ~ l/n> as shown in Fig. [3J It 
is important to emphasize the difference between finding 
scaling laws for averages and analyzing the worst case. 
From the point of view of characterizing a phase tran- 
sition, averages over Exact Cover instances follow quite 
well defined laws. The worst case is harder to discuss 
as no systematic search of it can be done. The worst 
case results can only be considered as consistent with 
the polynomial vanishing of the energy gap. It is worth 
noticing that the worst case, defined as the instance with 
a smaller minimum gap, brings also the higher entangle- 
ment as the system is passing closer to the phase tran- 
sition. Moreover, the minimum gap takes place at the 
same place where entanglement peaks. This phenomenon 
is illustrated in Fig^where s c appears to converge to the 
same value when n — > oo from above for the minimum 
gap and from below for the maximum entanglement. 

The scaling of entanglement entropy as the system 
passes near a quantum phase transition point is asym- 
metric, as is seen in Fig. ^ The growth of entangle- 
ment is slower in the beginning of the evolution and fits 
remarkably well a curve of the type E ~ log | log(s c — 
s)\, whereas the falling down of the peak is better 
parametrized by a power law E ~ |(s — s c)\~ a with 
a ~ 2.3. Both scaling behaviors improve as n becomes 
larger. 

Grover's algorithm ^l| does not produce exponential 
speed-up. It is then arguable that entanglement, despite 
being necessary for having some computational speed- 
up, should not play a relevant role in this case. It is 
indeed possible to analytically address this question. Let 
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Figure 3: Minimum gap versus the inverse size of the system, 
both for the worst case and for the average over all instances. 
Error bars give 95 per cent confidence intervals. The behavior 
appears to be linear. 
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Figure 4: Mean critical point s c for the minimum gap and for 
the maximum entropy. Error bars give 95 per cent confidence 
intervals. Both sets of points tend to approach as the size of 
the system is increased. 



us cast Grover's algorithm into the adiabatic evolution 
of the Hamiltonian 

H(s) - (1 - s)(I - \s)( s \) + s(I - \x )(xo\) , (8) 



where 



\x), n is the number of qubits, 



and \xq) is the marked state 0,0. Other alternative 
definitions of the Hamiltonian used in Grover's adiabatic 
algorithm (such as Hamiltonians cxplicitely defined in 
terms of qubits) do not lead to significantly different con- 
clusions from the ones presented in this work. The com- 
putation takes the quantum state from an equal super- 
position of all the possible computational states directly 
to the state \xo), as long as the evolution remains adia- 
batic. The time the algorithm takes to succeed depends 
dramatically on how we choose the parametrization of s 



in terms of time. The discussion on entanglement can be 
reduced to analyze its dependence on s since the explicit 
dependence on time and its consequences (see 0, for 
further information about this topic) are of no relevance. 

It is straightforward to check that the Hamiltonian © 
has its minimum gap between the ground and first ex- 
cited states at s = 0.5, which goes to zero exponentially 
fast as the number of qubits in the system is increased. 
Therefore, this Hamiltonian seems to undergo a quantum 
phase transition in the limit of infinite size as s = 0.5 be- 
comes a non-analytical point (for more on the Grover 
problem as a quantum phase transition, see We 
consequently expect quantum correlations to be maxi- 
mum for this value of s. We present without proof [l5| 
the result of the exact analytical calculation which shows 
that, for any equally sized bipartition, the entanglement 
entropy scales as 



E(s = 0.5, n -> oo) = 1 



2 o-n/2 

In 2 



(9) 



so the entropy tends to 1 for s = 0.5 as an square root 
in the exponential of the size of the system, which is the 
typical factor in Grover's quantum algorithm. Entangle- 
ment is also bounded for all possible bi-partitions and no 
exponential speed up is present. We must remark that 
our analysis is based on the study of the quantum state 
between succesive calls to the quantum oracle. Entan- 
glement in the quantum register might become very high 
during the application of the unitary black box, which 
will in turn depend on the specific searching problem we 
wish to solve. This very general situation can not be 
adressed in detail as it depends on the realization of the 
black box. As long as we restrict ourselves to the sit- 
uation between calls of the unitary oracle, we see that 
entanglement is a bounded quantity (which is not nec- 
essarily the case along the particular implementation of 
the quantum black box). 

The main theoretical challenge in quantum computa- 
tion theory remains quantum algorithm design. It has 
been observed that majorization theory seems to play 
an important role in the efficiency of quantum algo- 
rithms 0, Nevertheless, a relevant element for 
quantum computational speedup seems to be entangle- 
ment JfjUl^lirJ |2(j. Our results build on previous work 
EJ Sill and suggest that entanglement grows expo- 
nentially (as measured by x, the maximum rank of the 
reduced density matrices obtained over all possible bi- 
partitions) at a universal point along the adiabatic quan- 
tum evolution for the Exact Cover problem, which limits 
the possibility of an efficient classical simulation. At this 
point, the system comes close to a quantum phase transi- 
tion and entanglement obeys a scaling law that needs fur- 
ther investigation. For instance, it has also been proved 
that entanglement in Shor's factoting algorithm diverges 
exponentially fast in the number of qubits, which makes 
this algorithm difficult to simulate classically as well . 
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Scaling of entanglement seems to be further related to 
the effective connectivity of the system. Grover's prob- 
lem reduces to a two state problem and entanglement 
is bounded. The quantum Ising, XX and Heisenberg 
spin chains show logarithmic scaling. Higher dimensional 
spin models obey faster scaling laws. The maximum con- 
nectivity corresponds to non-local interactions (as those 
present in the adiabatic evolution algorithm for the 3- 
SAT problem) and entropy does approach its maximum 
scaling. Computationally hard problems are thus associ- 
ated to quantum systems that present phase transitions 
where entropy comes close to its maximum possible scal- 
ing. 
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